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Abstract 

We derive factorization identities for a class of preemptive-resume queueing systems, with 
batch arrivals and catastrophes that, whenever they occur, eliminate multiple customers present 
in the system. These processes are quite general, as they can be used to approximate Levy 
processes, diffusion processes, and certain types of growth-collapse processes; thus, all of the 
' processes mentioned above also satisfy similar factorization identities. In the Levy case, our 
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identities simplify to both the well-known Wiener-Hopf factorization, and another interesting 
factorization of reflected Levy processes starting at an arbitrary initial state. We also show how 
the ideeis can be used to derive transforms for some well-known state-dependent /inhomogeneous 
birth-death processes and diffusion processes. 
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■ 1 Introduction 

O 

^■f-^ . The Wiener-Hopf factorization is a classical result in both the theory of random walks and the 

theory of Levy processes. For a Levy process X, the factorization allows us to write the position 
of X at an independent exponential time e^, i.e. X(eq), as the sum of two independent random 
variables: info<s<e, X{s) and X{eq) — info<s<e, X{s), with the latter random variable representing 
. the reflection of X at a random time Cq. In principle, the distribution of the reflected process at 

time Cq can be derived if and only if the distribution of the infimum of X over [0, e^] is known as 
well. 

We show that a similar type of property is also found in processes that may not necessarily be 
expressible as a reflection of a simpler process. To do this, we introduce the Preemptive-Resume 
Production system, or PRP system, and we show that it satisfies a factorization identity. Technically, 
for an arbitrary PRP system the identity is not a true factorization, but it is in some cases: when X is 
a Levy process, for instance, our factorization identity is equivalent to the Wiener-Hopf factorization. 
The notion of a PRP system may appear at first to be somewhat contrived, but this is not the case: 
such systems can be used to approximate many types of important processes found in the probability 
literature, such as Levy processes, diffusion processes, and even Markovian growth-collapse models. 

Our factorization results also provide insight into the time-dependent behavior of a number 
of important birth-death processes, with birth/death rates that may depend on the state of the 
system. For instance, our Wiener-Hopf identity shows how the probability mass function of the 
M/M/s queue- length at an independent exponential time Cq can be expressed entirely in terms of 
quantities from a M/M/1 queue and a M/M/oo queue. Similarly, a M/M/s/K queue (assuming 
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s < K, otherwise trivial) can be expressed in terms of a M/M/oo queue and a M/M/1/{K — s) 
queue, and a similar observation may be made for Markovian queues with reneging. In particular, 
the pmf for the M/M/s/K queue can be quickly derived from the solutions to the M/M/oo queue 
and the M/M/1/{K — s) queue, without having to make use of the Kolmogorov forward equations 
corresponding to the M/M/s/K queue. Similar expressions can also be derived for diffusions that 
can be expressed as limits of birth-death processes. 

Readers wondering why we are interested in studying the distribution of X{eq) should note that 
P{X{eq) = k) can be expressed as q times the Laplace transform of the function P{X{t) ~ k) 
evaluated at q, where g is a positive real number. Hence, having knowledge of X(eq) yields insight 
into the behavior of X{t), for each t > 0. Even though we restrict ourselves to the case where q is 
real and positive, it is possible to derive similar transform expressions for the function P{X{t) ~ k) 
at complex numbers with positive real part: readers will find explanations of how to make such 
extensions at various places throughout the paper, whenever they are needed. 

The factorization results we present here seem to be somewhat related to those found in Millar 
[29] . The main result of [29] establishes that for a Markov process X satisfying suitable regularity 
conditions, the distribution of the path of X from the time at which a functional of it attains a 
minimum is independent of the behavior of X before having attained this minimum. Contrary to 
[29j . our factorization results are valid for processes that are not necessarily Markovian, and our 
results also show how various transforms associated with some processes can be decomposed into 
computable transforms associated with other types of simpler stochastic processes, as previously 
mentioned. 

2 Model Description 

We now define what we refer to as a Preemptive-Resume Production system, or PRP system. At 
time zero there are a countably infinite number of customers present, which are labeled no,no — 
1, no — 2, no — 3, . . .. The system then begins to process the work of the customer that possesses the 
highest label, or number, which at time zero is customer no. The server processes jobs in accordance 
to the Last-Come-First-Served Preemptive-Resume discipline. All customers possess a random, 
generally distributed amount of work, and the amount of work possessed by a given customer is 
independent of the amounts of work of all other customers that will visit, or have visited the system. 
We are interested in studying the process Q {Q{t); t > 0}, where Q{t) represents the label of the 
customer being served by the server at time t: for example, Q{0) — no. 

There arc two sets of Poisson processes governing arrivals to the production system. The first set 
governs single arrivals to the system, and consists of an independent collection of Poisson processes 
{Aojljgz, where Aqj has rate Xqj. At an arbitrary time t, when Q{t—) = j, we say that Aqj is 
active: in other words, if a point of Aqj occurs at time t while Q{t—) = j, then Q{t) = j + 1, and 
the new arrival is immediately given label j + 1. Otherwise, the point of Ao.j occurring at time t is 
ignored if Q{t—) = k j, so no new customer arrives to the system at that time. Once the server 
finishes with the customer having label j + 1, it begins serving customer j, returning to where it left 
off before previously departing. 

The second set of Poisson processes govern batch arrivals of customers to the system (we allow 
batches to be of size one). This second set consists of an independent collection of Poisson processes 
{Aij_k}j,k<£Z, where Ai^j^k has rate Xi_jP{Zi_j = k — j). Again, while Q{t—) = j, we say that the 
subcollcction {Aij^k}kei- is active, so a point of at time t pushes Q from level j to level fc, the 

k — j customers in the batch are instantaneously assigned labels j + I, j + 2, . . . , fc, and the server 
immediately begins processing customer k. Here Zi_j is a generic random variable representing the 
jump size of the Q process from level j: we allow the distribution of these jumps to depend on the 
current level. 

We further assume that catastrophes occur according to a modulated Poisson process D := 
{D{t);t > 0}, with rate (5g(t_). At the time of a catastrophe, a random number of customers are 
removed from the system: in particular, if Q{t—) = n, and a catastrophe occurs at time i, which 
eliminates k customers, then customers n, n — 1, n — 2, . . . , n — fc + 1 are immediately removed from 
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the system, and at time t the server begins to process the remaining amount of work possessed by 
customer n — k, and so Q{t) = n — k. We assume that the distribution function of the number of 
removals at time t depends on Q{t—), so that the downward jump distribution of the process may 
depend on the level of the process, immediately before a jump. 

Readers may wonder why we chose to use an infinite collection of independent Poisson processes 
to govern arrivals to our queueing system, while not modeling catastrophes in the same manner. The 
answer lies in the proof of our main result, as modeling the arrival processes in this way allows us to 
derive a linear system of equations in a most efficient manner. Indeed, catastrophes can be modeled 
in the same way, but these will not play as important a role in our proofs. Our use of collections 
of Poisson processes to model the arrival process was inspired by Chapter 9 of Bremaud [12] , who 
makes use of such a framework when constructing continuous-time Markov chains. Readers wishing 
to rigorously construct our PRP systems in the same manner can follow the procedure given there, 
by expanding the state space of the PRP system to include the residual service time of each customer 
in the system, thus making it a stochastic recursive system, and Markovian: readers should note 
that customers in the system possess generally distributed amounts of work, meaning {Q{t); t > 0} 
is not a Markov process unless the state space is expanded to include the residual service times. 

Later we will use these processes to approximate Levy processes: arrivals from the {^ojlj 
collection and service completions of the server will be used to construct Brownian motion, while 
the batch arrivals and catastrophe processes will be used to construct Compound Poisson processes. 

Finally, we also consider a 'reflected' PRP system {Qi{t);t > 0}, where I is a flxed integer. This 
system behaves in a similar manner as Q, with the following exception: whenever Qi is in a state i, 
and a catastrophe occurs which, in the original system, would place Q at a level at or lower than /, 
Qi instead makes a transition from state i to state I. When Qi is at level /, the server stops working 
until the next arrival: hence, customer I is in the system for all time. Finally, upward jumps of Qi 
behave the same as upward jumps of Q. We refer to Qi as a reflected PRP system with reflection 
at level 



3 Main Results 

Our main result establishes that the process {Q{t);t > 0} from the PRP system satisfles a factor- 
ization identity, which we now give. 

Theorem 3.1 Let Cq he an exponential random variable with rate q > 0, independent of Q. For 
any two integers fc, I, where k > and I < uq = Q{0), 

PiQi{eg) = k + l\Qi{0)^l) = P{Q{eg) = k + l\ inf Q{u) ^ I) 

= P{Q{eq) - inf Q{u) = k \ inf Q{u) = I). 



Proof To help readers understand the proof, we break it up into three steps. 

Step 1 We begin by presenting the following identity, which is satisfied by the sample paths of our 
PRP system: for each t > 0, we see that for any two integers k, I with k > 1, I < no = (5(0), 

l{Q{t) >k + l, inf Q{u) = I) 

0<u<t 

t 

l{Q{s-) = k-l + l, inf Q(M)=/)l(inf Qis) > k + l)Ao k^i+i{ds) 

0<u<s ue[s,t] 
k—l oo -.t 

+ yy / l{Q{s-)=j + l, inf Q{u)=l)l{mi Q{u) > k + l)Ai,j+i^,n+i{ds). (1) 

The identity (HI) says that, in order that Q{t) > k + I, exactly one of two things must happen: 
if the infimum of the process over [0,t] is I, either (i) there exists a time point s < t such that 
Q{s—) = k—l + l, Q{s) = k + l (due to the arrival of a customer from Aq at time s), and the process 
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stays at or above level k + I in [s,t\, giving the first term, or (ii) there exists a time point s < t 
such that, due to a batch of customers arriving at time s (which is contributed by ^i), the process 
crosses level k + 1, reaching some level at or above fc + Z at time s, and stays at or above k + 1 during 
[s,t\, giving the second term. 

After taking expected values of both sides of ([ij , we get 

P{Q{t) >k + l, inf Q{u) = /) 

0<lL<t 



Jo 

k—1 oo 



l(Q(s-) = fc-l + /, inf Q{u)^l)l{mi Q{u) > k + l)A„,k-i+i{ds) 

0<u<s u<£[s,t] 



"■ — ^ r nt 

+ YYE l(g(s-)=j + /, inf Q{u)^l)l{mi Q(u) > fe + 0^ij+;.m+!(ds) 

,, Uo 0<u<s ue[s,t] 



j—Q m—k 



(2) 



We can use the Campbell-Mecke formula to evaluate the expected values found on the right-hand 
side of Equation ([2]) . Notice first that 



E 



l{Q{s-)^k-l + l, inf Q{u)^l)l{ inf Q{u) > k + l)Ao k-i+i{ds) 

0<u<s ui£[s,t] 



Xa.k-i+i Vs{Q{s-) = k-l + l, inf Q{u)^l, inf Q{u)>k + l)ds, 

0<u<s ue[s.t] 



where V represents the Palm kernel induced by Ao.fc-i+;. Furthermore, since the server processes 
work in a preemptive-resume manner, we can also use the Campbell-Mecke formula to establish that 

T'.C inf Q{u)>k + l,Q{s-) = k-l + l, inf Q{u) ^ I) 

ue[s,t] o<«<s 

= P{Tk+i,k+i>t^ s)Ts{Q{s-)^k-l + l, inf Q{u) = l) 

0<u<s 

where Tkj is the amount of time it takes the PRP system to go below state j, starting from state 
k, j < k, where all customers labeled j, j + 1, . . . , k have not yet received any attention from the 
server. Moreover, if we let {J-j; < > 0} represent the minimal filtration induced by Q and our arrival 
and catastrophe processes, we see that the event {Q{s—) = k — l + l, info<„<s Q{u) = 1} ^ J^s~, and 
so Proposition I A. II in the Appendix yields 

rs{Q{s-) = k-l + l, inf Q{u) = l)=P{Q{s)=k-l + l, inf Q(u) = 0- 

0<^i<s 0<u<s 

An analogous argument can be used to evaluate the second type of expectation found in ([2]) . Plugging 
these expressions into ([2]) gives 



0<u<t 



P{Q{t)>k + l, inf Q(u) = = Ao,fe-i+i / P{Q{s)=k-l + l, inf Q{u) ^ l)P{Tk+Lk+i > t - s)ds 



0<u<s 



k—1 oo 



+ V V Xl.j+lP{Zi,,+i = 771 - j) / P{Tm+l,k+l > t - s)P{Q{s-) ^j+l, inf Q{u) = l)ds. (3) 
' If, 0<u<s 

j=Om=k - - 

After integrating both sides of ([3]) with respect to an exponential density with rate g > 0, we get 

(1 - (pk+Lk+i{q)) 



P{Q{eg)>k + l, inf Q{u) = I) = Xa,k-i+i- 

k— 1 oo 



.1^ (j.^(_q 



q 



-P(Q(e,) = fc - 1 + /, inf Q(it) = /) 



0<«<e, 



j—0 m—k 



; (1 - 0: 



m+Lk+l 



(q)) 



q 



PiQieq)=J+l, inf Q{u)^l) 

0<u<en 
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where (f)rn+i,k+i represents the Laplace- Stieltjes transform of T.m+i.fc+;(0) (with (5(0) ~ m + I). Di- 
viding by P(info<u<e, Q{u) = I) finally yields 

P{Qieg)>k + l\ inf Q(ii) = /) = Ao,fc-i+i-^i^^^±^^^±^P(Q(eg) = fe-l + M inf Q{u) ^ I) 

0<u<eq q 0<u<eq 

+ E E Ai.,+^P(Zi,,+, =m-j) ^^~'^'"+"+'^'^^^ P(Q(e,) =j+ / I inf Q(u) = 0- (4) 

Q 0<u<e„ 
j=0 m=k ^ - - 9 

Step 2 We now show that the system of equations Q has a unique solution. Notice that for a fixed 
integer these equations can be iteratively solved, since 



y P{Q{e,) = k + I \ inf Q{u) = l) = l. 

k=0 - - 1 



Indeed, notice that 



1 - P{Q{eq) = l\ inf Q{u) = l) = P{Q{eg) >l + l\ inf Q(s) = /) 

0<u<eq 0<u<eq 

^ xJl^^I±h!±lMp^Q^e,)=l\ inf Qiu)^l) 

q 0<M<e, 

•■^ — ' q 0<«<e„ 

m— 1 

which allows us to determine P{Q{eq) = I \ info<ti<e, Q(u) = 0: other probabilities can 

be determined in a similar, iterative manner. Hence, there is a unique probability measure on the 
integers that satisfies these equations. 

Step 3 By precisely the same arguments, we see that the Qi process satisfies the same system of 
equations. Indeed, when Qi{0) = I, 

PiQiie,) >k + l) = Ao,fe+i_ii^^^±^^^±^P(g,(eJ = k-l + l) 

q 

+ E E Au+,P(^u+, = - - j)i^^i^i±i^p(gKe,) = J + I). 

j=0 m=k 

Thus, we see that 

P{Qi{eq) = k + l\Qi{{))^l)^P{Q{eq)=k + l\ inf Q{s) ^ I) 

0<s<eq 

completing the proof. <^ 

Remark It is worth noting, from the point of view of numerical transform inversion [4], that a 
similar result can be derived when we consider complex-valued q, i.e. expressions of the form 

P{Q{t) ^k + l, inf Q{s) = l)qe~'i*dt 

0<s<t 

for complex q with positive real part, i.e. those q satisfying 3fi(g) > 0, as opposed to P{Q{eq) = 
fc -f /, info<s<e5 Q{s) — I) for real q > 0. First note that for g = x -I- ly satisfying 3fi(q) = x > 0, with 
Cx being exponential with rate x, independent of Q, 

P{Q{t)^k + l, inf Q{s) ^ l)qe~'^*dt = / P{Q{t)^k + l, inf Q{s) = l){x + iy)e-'y*e-'-'*dt 

0<s<i /n 0<s<t 



^^i^P[l(Q(e,) = fc + /, inf Q(s) = Oe^'^^^ 

0<s<ej: 
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Using this observation, wc can mimic the proof of Theorem 13.11 in a straightforward manner to 
determine that 

E[l{Q{e,) = k + l)e-'y^- \ inf Q{s) = I, Q{0) = no] 

0<s<e-^ 

E[e-'y'- l(info<,<e, Q{s) = I) \ QjO) = no]E[l{Qi{e,) = k + l)e-'y^- \ Q;(0) = l]x + iy 
P(info<.<e, Q(s) = I I g(0) = 7io) X 

which contains quantities that are given in terms of either the reflection Qi reflected at I, or hitting- 
time transforms associated with the original process Q. To see why only these types of transforms 
need to be computed, note that letting t; = vn.i{t > : < /} yields 

E[e~'y^-1{ inf Q{u) = /) | g(0) = no] = Ele-'^^^li inf Q{u) < I) \ Q(0) ^ no] 

- i;[e-*^'=n( inf Q{u)<l-l)\Q{0)=no] 

0<u<ex 

= E\e-'y^^\{Ti<e^)\Q{^) = no] 

- i;[e-*^'=n(Tz_i < e,) I Q(0) = no] 

= ^—E\e-'y^'\(Ti < e.) I Q(0) - no] 

- ^-E[e-'y^'-n{Ti^, < e,) I Q(0) = no] 

X + ly 

= I 0(0) = no] - i?[e-«-'- | Q{0) = no]] 

X + ly 



E[l{Q{e,) + l)e~'y^- \ inf Q{s) = I, Q{Q) - no] 



This gives 

[g[e-^"' I Q(0) ^ no] ~ gje-g^'-i | Q(0) ^ no]] 

[E[e-^^^ I g(0) = no] - E[e~--^-^ \ Q(0) = no]] 

implying 



E[l{Qi{e^) = k + l)e-'y-^ I Oi(0) = Z] 



P{Q{t) = k + l, inf Q{s)^l\Q{0)=no)qe~'^^dt 

0<s<t 

[E[eri^' I 0(0) = no] - E[e'''^'-^ \ 0(0) = no]] H P{Qi{t) = k + l\ 0/(0) ^ l)qe~^'dt 



which is clearly the complex analogue of the formula given in Theorem 13.11 All other types of 
transforms that we will need can be computed in a similar manner, for complex q. 

We now show that the reflected process {Qo{t);t > 0} exhibits a similar type of factorization 
identity. 

Theorem 3.2 Suppose Q is a PRP system with 0(0) = no, and let Oo be the reflected version of 
at level zero, with Oo(0) = no. Then for each integer I > 0, and each integer k > 1, 

P{Q{eq)- inf Q{u)^k\ inf Q{u) = I) = P{Qo{eq) - inf Qo{u) ^ k \ inf Qo{u) = I). 

0<u<eq 0<u<eq 0<u<.eq 0<(i<eq 

Proof Notice that a sample-path identity that is completely analogous to ([l]) can be established 
for Oo- for each I > 0, k > 1, 

l(Oo(i) >k + l, inf Oo(m) = 

0<u<t 

= / l(Oo(s-) = fc-l + Z, inf Oo(") = 0l(inf Qo{u) = k + l)Ao,k-i+i{ds) 

Jq 0<u<s ue[s,t] 

k—1 oo pt 

VV / l{Qo{s-)=j + 1, inf Qo{u) = l)l{ inf Qo{u) > k + l)Ai,j+i,^+i{ds). (5) 

I. Jo 0<u<s u<£[s,t] 



j—0 m—k ' 
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Applying the same steps found in Step 1 of the proof of Theorem 13.11 yields 

P{Qo{eg) > k + I \ inf Qoiu) ^ I) ^ Xo.k-i+i ~ ^'^^^ P{Qo{eg) ^k-l + l\ inf Qo{u) ^ I) 

0<u<eg q 0<u<eq 

+ J2 E Ai,,-+^P(Zi.,-+, ^m~j) ^^~'^'"+'-"+'^'^^^ P(go(e,) inf Qoiu) ^ I). (6) 

j—0 ni—k 

For our fixed I, we notice that the equations that form system ^ are the same as the equations 
found in Hence, by the uniqueness result proven in Step 2 of Theorem 13. II we have 

P{Q{eq)>k + l\ inf Q{u) = I) = P{Qo{eg) > k + I \ inf Qo{u) ^ I) 

which completes the proof. <C> 

Two interesting factorization results can be derived, when the batch and catastrophe sizes of 
both Q and Qq have distributions that are state-independent. Clearly, in this case we see that for 
each fc > and I, P{Qi{eq) = k + l \ 0/(0) = = P{Qo{eq) = k \ Qo{0) = 0), and since Qo is the 
reflection of Q at level 0, we also find that 

<9o(eq) = Qieq) - inf Q{u) 

which follows since customers are processed in a Last-Come-First-Served Preemptive- Resume man- 
ner. Hence, Theorem 13.11 vields for each k > 0, I < Q = Q{0), 

P(0(eg)- inf Q{u) ^ k) ^ P{Q{eq) - inf Qiu) = k \ inf Q{u) ^ I). 

0<u<eg 0<u<eg 0<u<eg 

In other words, the following corollary holds. 

Corollary 3.1 Suppose that {Q(t);t > 0} represents a PRP system, with state-independent jumps, 
and let Cq be an exponential random variable with rate q > 0, independent of Q. Then for each 

Here is the expectation corresponding to Px, where Px is a probability measure under the 
condition that our process starts at level x. This notation will be used in many places throughout 
the rest of the paper. 

This factorization has been well-known for Levy processes since the late 60 's, due to Pcrcheskii 
and Rogozin |31| , and the first probabilistic proof of this result was given in Greenwood and Pitman 

We can also conclude from Theorem 13.21 that for / > 0, when Qo{0) = Q{0) ~ uq, 

P{Qo{eq)- inf Qo{u)^k\ inf Qo{u) = I) = P{Q{eq) - inf Q{u) = k) 

= P{Qo{eq)- inf goH = fc) 

where the second equality follows from the simple fact that the reflection of Qo at its infimum is equal 
in distribution to the reflection of Q at its infimum. Hence, we see that (5o(e<j) — info<u<e, Qaiu) is 
actually independent of info<ti<e, Qa{u), which gives us another interesting corollary. 

Corollary 3.2 Suppose that {QQ{t);t > 0} is a reflected version of our PRP system, reflected at 0. 
Then for each cj G M, and each integer uq > 0, 
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Such a factorization result is useful when studying reflected processes starting in an arbitrary 
initial state. Corollarv l3.1l shows that, since info<u<e^ Q{u) is independent of Q{eq) — mio<u<eg Q{u), 
the transforms oiQ{eq) and info<u<e<j Qiu) can be used to derive the transform of (5(eq)—info<M<eg Q{u), 
which represents the distribution of the reflected process, starting in level zero. Theorem 13.21 can 
then be used to find the distribution of the reflected process, starting in any initial state, since it is 
clearly equal in distribution to a convolution of the reflected PRP system Qg starting in level zero, 
and a truncated version of info<M<e Q{u). 

We are now ready to see how the Wiener-Hopf factorization for Levy processes follows as a 
consequence of our factorization identities for PRP systems, whose arrival rates, service rates, and 
jump distributions do not depend on the level of the process. 

3.1 The Wiener-Hopf factorization 

We begin with establishing the well-known version of the Wiener-Hopf factorization, for Levy pro- 
cesses. 

Theorem 3.3 Suppose X is a Levy process, and let Cq he an exponential random variable, indepen- 
dent of X , with rate q > 0. Then info<s<e, X{s) and X{eq) — info<s<eg X{s) are independent. 

Proof Suppose first that X is a Levy process that consists of only a Brownian component and a 
compound Poisson component. In this case, there exists a sequence of PRP systems {X„}„>i, such 
that Xn converges uniformly on compact sets to X: in fact, each X„ process is also a Levy process. 
We omit the details on constructing the {X„}„ sequence, as they are somewhat standard: interested 
readers can also find them in a previous online version |19| of the paper. 

From Corollary we see that the Wiener-Hopf factorization is valid for each PRP system with 
state-independent jumps. Applying the Levy continuity theorem yields, for each (101,102) G K^, 

info<s<e, X(s)+a;2(X(e,)-info<s<e, X{s)))^ ^ £'[e'("i info<s<e5 ^„ (s)+tJ2 (e<,)-info<^<eg X„(s)))j 

n— f 00 

= iJ[e*"^ iiifo<s<s5 -^(s)]^Jgjt^2(^(e,)-info<s<e, X{s))j 

proving independence. To derive this result for an arbitrary Levy process, use this result in con- 
junction with the proof of the Levy-Ito decomposition: again, finer details of this procedure can be 
found in [H]. 

Our idea of proving a factorization result for a special type of process, then taking limits is similar 
to the older approaches of proving the Wiener-Hopf factorization, along with related results: see for 
instance Percheskii and Rogozin [3T], along with Gusak and Korolyuk [21]. Our approach differs 
in the fact that we use a discrete state space in continuous time: this allows us to state a simple 
sample-path identity, from which we derive a linear system of equations that has a unique solution. 
Moreover, our limiting argument makes use of classical heavy-traffic results from queueing theory. 
Readers interested in learning more about classical approaches towards proving the Wiener-Hopf 
factorization are referred to the recent paper of Kuznetsov |27| . 

3.2 An analogous factorization for the reflection 

We now show how to use Corollary 13.21 to deduce an analogous factorization for reflected Levy 
processes, with an arbitrary initial state. 

Theorem 3.4 Suppose X represents a Levy process, and let Cq be an exponential random variable 
with rate q > 0, independent of X. Moreover, let R := {R(t);t > 0} represent the reflection of X, 
with a reflected barrier at state zero. Then, assuming X{Qi) ~ x > 0, 
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Proof The proof of this rcsuh is completely analogous to the proof of Theorem 13.31 First, we use 
Corollary 13.21 to establish that it holds for a Levy process X that consists of only a Brownian and 
compound Poisson part. The general statement then again follows as before, from the proof of the 
Levy-Ito decomposition. 

Theorem 13.41 can also be derived directly from the Wiener- Hopf factorization. Here X{0) — x, 
and for each t > 

inf min(X(s),0) 

0<s<t 

and so 

Rit)- M R{s) = X{t) - M mmiX{s),0) - mi ({Xis) - mi min(XH,0) 

Let To = inf{< > : X{t) = 0}. If tq > t, then 

R{t)- mi R{s) = X{i)- mi X{s) 

since min(X(s), 0) = for < s < tq. Next, if tq < t, we also see that 



R{t) - inf ,i?(.s) = X{t) - inf X(,s) - inf X{s) - mi X{u) = X{t) - inf X(,s) 

0<S<t 0<S<t TO<S<t \ TQ<U<S I 0<S<t 

since miTg<s<t {X{s) — mir„<u<s X (u)) > 0, and X{to) — infT-,-|<„<T-,| X{u) = 0. Moreover, for each 
t > 

inf R(t) = max( inf X(s),0). 

0<s<t 0<s<t 

Thus, for an exponential random variable with parameter q > 0, independent of X, we have 

/>oo 

^^Jgj"(-R(e5)-info<s<c, i?.(s))git^info<s<eg i?.(s)j _ / ^g4i<j(-R(t)-info<3<t fl(s))giw info<s<t i?,(s)j^g-(jt^^ 







= £;^[g'"(-''^('^<!)~'"fo<s<=g ('*))gi'^rnax(0,info<,<e, X(s))j 

= i5^[g'"(-^('=<!)-info<s<eg X(s))j^_^jgiLjmax(0,info<s<e, X{s))j 



where the last step follows from the Wiener-Hopf factorization, i.e. Theorem [ 

Theorem [nH] does not seem to be explicitly known, however direct computations of E^le^'^^'''^''^] 
have appeared in various places: see e.g. Theorem 9.1 of Abate and Whitt [T], Theorem 2.1 of 
Abate and Whitt [3|, Bingham [10|, Bekker et al. [9], and Chapter 9, Theorem 3.10 of Asmussen [6], 
where all of these references address the factorization in the case where X is spectrally positive, i.e. 
X has only positive jumps. Theorem 13.41 is also implicitly stated in Example 3 of Palmowski and 
Vlasiou [30], in terms of the steady-state distribution of a reflected Levy process that experiences 
catastrophes at times forming a homogeneous Poisson process. Their result, like previous references, 
considers only the spectrally positive case, but their arguments can also be used to establish Theorem 
3.4 as well. Other results similar to Theorem 13.41 can also be found in the recent work of Debicki et 
al. [ini, and in Kella and Mandjes |26) . 



4 Applications to birth-death processses, and diffusions 

We now apply our factorization identities, i.e. Theorems 13.11 and 13.21 towards the study of birth- 
death processes, which form another interesting subclass of PRP systems. It will also be possible to 
apply our identity towards the study of diffusion processes as well, as these are often weak limits of 
birth-death processes. 

Readers should note that the transforms derived below can also be modified so that the domain 
is complex- valued, as we noted in the remark following Theorem 13. II above. 
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4.1 Birth-death processes 

Suppose that Q := {Q{t)\t > 0} represents a birth-death process on the mtegers, with birth rates 
{Anjngz and death rates {/inlnsz- Let represent an exponential random variable with rate q > 0, 
independent of Q. Throughout we assume that Q is ergodic, and we let tt represent its stationary 
distribution. Our object of study is now the probability mass function of Q{eq). We remind readers 
that Q can easily be related to a PRP system: units arrive according to a collection of independent 
Poisson processes {^ojljez where Aqj has rate Xj, each customer brings to the system a unit 
exponential amount of work, and the server processes work at a rate /i„ whenever the system is in 
state n, for n € Z. 

By Corollary 4.1.1 of Abate and Whitt [3], we see that for each n G Z, 

Po{Q{eq) = n) = 



where P„ is meant to represent a conditional probability, given Q(0) = n. This expression also holds 
in the absence of ergodicity, and also for complex q when Po((3(eq) = n) is interpreted as a Laplace 
transform, multiplied by q. 

However, suppose we would like to change the initial condition. While the same method will tell 
us that 

for an arbitrary no, we must be careful: how do wc know that £'„[e~*'^"o] is tractable? This is a 
very legitimate question, as there are many instances where £'„[e~^'^""] will be tractable for some 
choices of ng, but not for others. 

Thus, the key to computing these probabilities is to choose the appropriate reference point, i.e. 
the point found in the hitting-timc Laplacc-Sticltjes transforms given in the pmf of Q{eq). This is 
where our factorization identities become useful: they allow us to use whatever reference point we 
like, regardless of the initial value. 

We illustrate our approach by computing the pmf of the number of customers in an M /M / s 
queueing system at an independent exponential time Cq. The reader will see that our expressions 
will be given in terms of an M/M/1 model and an M/M/oo model, which are much simpler. 

4.1.1 The M/M/s queue 

Recall that the M/M/s queue is a birth-death process on {0, 1,2,.. .} with birth rates A„ = A, for 
71 > 0, and death rates ^„ = min{n, s}/x, for n > 1. A classical reference on the time-dependent 
behavior of the M/M/s queue is Saaty [35], which makes use of the approach found in Bailey [5]. 
Assume first that Q(0) = s. In this case, for each n > 0, 

P{Q{e)-n)- ""^"[^'^"] 

This is a nice expression: notice that if fc < s, Ek[e~'^'^^] is the Laplace-Stieltjes transform of the 
amount of time it takes an M/M/s queue to go from level k to level s, but this is the same as the 
Laplace-Stieltjes transform of the amount of time it takes to go from k to s in an M/M/oo queue, 
with arrival rate A and service rate ^. Similarly, for k > s, Ek[e~'^'^''] is just the LST of the amount 
of time it takes to go from level k to level s in an M/M/\ queue, with arrival rate A and service rate 
s[i. Hence, all of the terms in our expression for Ps{Q{eq) = fc) can theoretically be derived from 
two simpler models, the M/M/1 queue and the M/M/oo queue. 

For fc > s, we already have a closed-form expression for Ek[e~'^'^'']: letting -0(g) ~ E'g+i [6^^"^=] be 
the busy period of an M/M/1 queue with arrival rate A and service rate s/i, we see that 
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Wc now focus on the case where k < s. Letting {Qa//a//oo(0! ^ ^ 0} represent the queue-length 
process of an M /M/oo queue (including the customers in service), we use a classical argument found 
in Darling and Siegert [TS] to find that 

-Pfc (Qm/M/ oo (eg) = S) = Pk{QM/M/oc.{eq) = S,Ts < Cg) 

= PsiQM/M/M^s)Ek[e~^-^] 

giving 

p r -gr,i Pk{Q M /M /ocjeq) = s) 
Ps[QM/M/oo[eq) = s) 

To compute Pk{QM/M/oc,{^q) = s), we need to use the following known lemma. The /i = 1 case 
was observed in Flajolet and Guillemin |17| . but we repeat it here for convenience. 

Lemma 4.1 For a positive real number q, 



oo 

qe- 



where M is Kummer's function, i.e. 

with (a)o ~ 1, and for n > 1, (a)„ — (a)(ci + 1) ■ • • (a + n — 1). 
Proof Applying partial integration gives 

Jo 

After repeatedly applying partial integration and taking limits, we get the result. <^ 



Lemma 4.2 For each k < s, 

Proof This identity can be derived from the known fact that, at a fixed time i > 0, Q{t) is the 
convolution of a binomial random variable with parameters (fc, e~^*) and a Poisson random variable 
with parameter p(l — e"''*). The result then follows by integrating the pmf of Q(t), and applying 
Lemma l4Tl 

By making use of this lemma in equation ([8]), we arrive at the following result. 
Lemma 4.3 For each k < s, we see that 

E^^o E':::o'' o (i, g + , + ^ + 1, -p) 



Remark As discussed in the remark following Theorem 13.11 Lemmas 4.1, 4.2 and 4.3 can be 
modified so that q is allowed to take on complex values. 

Our next step is to use the Wiener-Hopf identity to compute probabilities of the form Pk{Q{eq) = 
7i), for arbitrary k,n>Q. Notice that we already have a nice expression for such a pmf, when k ~ s. 
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Case 1: k > s, n < s. Notice that 

PkiQisq) = n) = Pk{Q{eq) ^n,Ts < eg) 

= PkiQieg) - n I < eg)Ek[e-''^^] 
= P.(g(e,) = n)Sfc[e-?-=] 

showing, from our previous calculations, that this probability is tractable. Readers should again 
note that a similar argument can be made for complex q = x + iy satisfying x > 0. Here 

'-Ek[e~'y^-'^-^^+^^h{Q{e^ - + r,) = n)l(e, > r,)] 
-S,[e-'^^n(g(e,) = n)]Ek[e~'y-n{T, < e,)] 



X 


+ iy 




X 


X 


+ iy 




X 


X 


+ iy 




X 


X 


+ iy 


X 



Jo 



where the fourth equality holds by the strong Markov property. 

Case 2: k > s, n > s. This case is much more interesting, since it is possible for our process to go 
from k to n, without ever reaching level s in [0,6^]. Proceeding in the same manner as in Case 1 
yields 

Pk{Q{eq)=n) = Pk{Q{eq) ^ n,Ts <eq) + Pk{Q{eq) ^ n,Ts> Cq) 

min{n,A:} 

= Ps{Q{eq) = n)Ek[e-'^'']+ V Pfc(Q(e,) = n | inf Q{u) = l)Pk{ inf Q{u) 



0<u<eq 0<u<eq 
l—s-\-l 



These terms are computable: first note that 



Pk{ inf Q{u) = l) = Pkin < eq) - Pi^in.i < e,) 

0<.u<.eg 

= Ekie-"^'] - Ek[e-'^^^-'] 

and from Theorem 13.21 we find that conditional on info<ti<e, Q{u) — I, Q behaves as an M/M/1 
queue on [0, Cq] with arrival rate A and service rate sfi. Hence, 

\^/j{q)\ fXiiqy-' 



PiiQieq) = n\ inf Qiu) = I) = { 1 - 

0<u<eg \ S/i y \ S/i 

Case 3: < k < s, n > s. This case is analogous to Case 1: here 

PkiQieq) = ?i) = PsiQicq) = n)Ek[e''>^^^]. 

Now we can use Lemma 14.31 to express i?fc[e~'^^=] in terms of Kummer functions. 

Case 4: < k < s, n < s. As expected, this case is analogous to Case 2, but the expression here is 

more complicated than the other cases. Here 

Pk{Q{eq)^n) = Pk{Q{eq) ^ n,Ts < eq) + Pk{Q{eq) ^ n,Ts > eq) 

= P,{Q{eq) = n)Ek[e-''^^] + ^ Pfc(Q(e,)=n| sup Q{u) ^ l)Pk{ sup Q{u) 

l-maxikn\ 0<u<e, 0<u<e, 
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However, we again observe that 



Pk{ sup Q{u)^l)^Ek[e-'^-']-Ek[e~'^-'+'] 

0<u<-eg 

and conditional on supo<u<e^ Qi^) — I, we use TlLeorem l3.2l to deduce that Q behaves as an M/M/l/l 
queue on [0,6^], starting at level I. This yields 



Pk{Q{u) = n I sup Q{u) = I) 



implying that this final case is tractable as well, in that it can be expressed in terms of Kummer 
functions. 

There is an important lesson to be learned from our calculations of the pmf of Q{eq). Given 
a proper choice of initial point and reference point, our probability mass function of Q{Eq) can be 
expressed in terms of quantities related to three simpler models: the M/M/1 queue, the M/M/l/l 
queue, and the M /M/oo queue. Had we chosen another reference point different from s, our hitting- 
time transforms would have been much more diSicult to compute. 

4.1.2 The M/M/s/K queue 

Our factorization identities can also be used to derive the pmf of the M/M/s/K queue-length process 
at an independent exponential time e^, where s is the number of servers and K the system capacity. 
By choosing our reference point to be s, we mimic the procedure used in the M/M / s case to express 
the desired pmf in terms of two simpler models: the M /M / s/ s queue (which is expressible in terms 
of M/M/oo hitting-timc transforms), and the M/M/1/{K — s) queue. 

Note that the relevant hitting-time transforms for the M/M/1/{K — s) queue can be derived 
from the M/M/1 queue, since we can use the pmf of an M/M /I queue at an exponential time to 
derive the LST of the time it takes us to go from level ji to level j2 in an M/M/1 queue, when 
ji < j2 • Such a result can then be used to derive all of the corresponding hitting-time transforms 
for an M/M/1/{K - s) queue. 

4.1.3 Time-dependent moments 

It is possible to make use of the factorization identities to derive the moments of Q{eq) as well. 
To illustrate the main idea, we first suppose that {Q{t)]t > 0} represents an M/M/1 queue-length 
process, with arrival rate A and service rate ji. It has been shown in Abate and Whitt [2] that, for 
each < > 0, 

EMt) I Q(0) = 0] = -^P{Rr < t) 

where r represents the busy period of an M/M/1 queue, and Rr represents the residual busy period, 
i.e. for each t > 0, 

P{Rr >t) = -^ f P{t> x)dx. 
Letting be an exponential r.v. with rate q > 0, independent of Q, gives 

- P 

p 1 - E[e~i^] 



E[Q{eq)\Q{{)) = Q\ = 



1 - p qE[T] 
A(l - Eje-"^]) 

q 
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which imphes that the first moment of Q{eq) is tractable, assuming we start in state 0. 

Our factorization identities can now be used to compute the first moment of Q{eq), for any initial 
condition. Suppose that Q(0) — no >0. Then 

E[Q{eg) I Q(0) = no] = E[Q{eg) \ inf Q{s) ^ 0, Q{0) = no]P{ inf Q{s) ^ | Q{0) ^ m) 

O^S^Cq 0<S<eq 



+ E^[QK)L>f Q(s)-fc,Q(0) = no]P( inf Q(,s) = fc | g(0) 

fe=0 

= E[Q{eg) I Q(0) = 0]P{ inf Q{s) = | Q(0) = no) 

0<s<eq 

+ YiE[Q{eg)\QiO) =0]+k)P{ inf Qis) ^ k \ Q{0) = no) 
fc=o 

Wo 

= £;[Q(e,) I Q(0) = 0] + y kP{ inf Q(.s) = fc | g(0) = no) 

^ — ^ 0<s<eo 
A;^0 

= ^^i-^i^+i:fc^(gr-'=(i-^(<z)). 



fc=i 



The key step in this derivation is the second equality: if info<s<e ~ fc, then Theorem 13 . 2 1 tells 

us that Q(e<j) is equal in distribution to an M /M/1 queue on the states {fc, fc + 1, fc + 2, . . .} with 
arrival rate A and service rate fi. This result agrees with the result given in [3], and also in |20j . 
With a bit of patience, higher moments can also be computed through the use of this approach, but 
there are better ways to do this for the M/AI/1 model: see |2D] for details. 

An analogous procedure can be used to compute the moments of Q{eq), for more complicated 
processes. Suppose now that {Q{t)] t > 0} represents the queue-length process of an MjMj s queue, 
with arrival rate A and service rate /i, and s servers. While the transient moments of the AI /M / s 
queue have been studied in Marcellan and Perez |28j , the point here is to show how to construct the 
moments from simpler birth-death processes. 

The key to computing the moments of Q{eq) for an arbitrary initial condition is to first compute 
the moments, while assuming that (3(0) = s, since we will want to again use s as a reference point 
when we apply Theorem 13.21 Again, since Q is a reversible process, we can say that 



k 

k—s 



fc=0 k=s+l 



E[Q{eq)\Q{0) = .s]=7:o{q)Y.kEk[e-"i^^]f^+no{q)'^ E kEu[e-'i^^]{p/ s 

with 

■T^a{q) = 



.fe=0 k=s+l 



being the normalizing constant. There are a few observations here worth noting. First, notice that 
7ro((7)Efc^fc[e"^"lf[ = Ps{Q{eg)<s)E,[QM/M/s/s{eg)] 

where Qm/m/s/s represents an M/M/s/s loss model with arrival rate A, service rate p, and s servers, 
and this is a known expected value; see Abate and Whitt |5] for details. Second, we see that 

g oo g oo 

Mq)^ E kE,[e-''^^]{p/s)'^-'' = no{q)^ E (k - s)E,[e-'^-^]{p/ s)''-^ 

fc=s-|-l ' k=s+l 



+ Mq)^ E ^Ekie-o-^Kp/s)"-" 

k=s+l 

= Ps{Q{eq) > s)Eo[QM/M/iieq)] + sPs{Q{eq) > s)Po{Qm/m/M > 1) 
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where Qm/m/1 represents an M/M/1 queue with arrival rate A and service rate s/i. Thus, we 
conclude that E[Q{eq) \ (3(0) = s] is a quantity that can be computed. 

To get E[Q{eq) \ Q{Q) — i] for an arbitrary i > 0, we now invoke Theorem 13.21 Suppose first 
that i < s. Then 

s-l 

E[Q{eq) I Q(0) =i] ^ ^[^(^9) I s'^iP ^ J' = *]^( s'^iP Q(^) = ^ I = 

■ ■ 0<s<eg 0<s<eq 

+ i;[Q(e,) I Q(0) = s]P(r, < e,) 

and we observe from Theorem 13.21 that, conditional on supQ<^<g^ Qis) — J, Q{&q) behaves as an 
M/M I i / j queue on {0, 1, 2, . . . , j}, meaning 

E[Q{eq)\ sup Q{s) = j,Q{Q) =i] = E.j[QM/M/j/j{e.q)\- 

0<s<e, 

All of the other terms in the sum are, for similar reasons, also tractable. A similar argument can be 
used to derive E[Q{eq) \ (5(0) = i] for i > s; we omit the details. 

We also point out that a similar argument can be used to derive moment expressions for the 
M/M/s queue with exponential reneging, i.e. the M/M/s — M queue, which is the model studied 
in Garnett et al. [22] ■ Such moments would be decomposed into components from an M/M/s/s 
queue, and a M/M/1 — M queue, and the M/M/\ — M queue moments have recently been studied 
in [H]. 



4.2 Diffusion processes 

The factorization identities can also be used to establish similar expressions for diffusion processes. 
We illustrate how the procedure works by applying it to a classical reflected diffusion: regulated 
Brownian motion. 



4.2.1 Regulated Brownian motion 

Suppose that {B{t);t > 0} represents a Brownian motion, with drift /z = — 1 and volatility ct^ = 1. 
We are interested in understanding the time-dependent behavior of {R{t);t > 0}, where 

R{t)^B(t)- inf mm{B{u),Q) 

0<lL<t 

i.e. R is the one-sided reflection of B. Granted, since S is a Levy process, we can already use 
the Wiener-Hopf factorization to derive the Laplace-Stieltjes transform of R{eq). However, we will 
instead be interested in showing how our factorization identities can also be used to derive the 
probability density function of R{eq). 

To derive this pdf, we will need to know a bit about the distribution of the hitting times associated 
with a Brownian motion. Following the classical argument of applying the optional sampling theorem 
to the Wald martingale, we see that 



E,.[e- 



Moreover, R has a unique stationary distribution tt, where Tr{dx) = 2e''^^dx. 

We will now compute the density of R{eq), given i?(0) = xq: we denote this density at the point 
X as fji(^e )(x; xq). Again, we will need to break the calculation up into cases. Considering first the 
case where x > xq, we may use Theorem l3.2[ along with a weak- convergence argument to show that 

P,JR{eq)>x) - iJ^Je-9-o].£:^?^l£l!I!MM 

""^ Ey\e~'^'-"]7T(dy) , , . . 

^-^1 r-T-^rf^' inf RW) < z). 

Ey[e-i-'o]TT{dy) V«<e, ' - ' 



f Jx 

poo 
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Careful readers will note that this identity is valid for a large class of reflected diffusion processes 
(namely, those processes that are expressible as a scaling-limit of a sequence of birth-death processes), 
not just for regulated Brownian motion. Success in using this identity for a given diffusion depends 
on both the tractability of the hitting-time transforms, and the integrals containing them. 

For a; > 0, we can use our expressions for both the hitting-time LST and the stationary distri- 
bution to show that 



/•oo 

J X 



]TT{dy) 



g-(-i+vT+2^)y2e-2y^y 
2 



l + y/TT2q 



Also, for < z < xo. 



inf R{u)<z) 



-(-l + VTT2?)(xo-z) 



so for positive z, we find that the density of info<ti<e, R{u) is just 



dP{ inf i?(w)<z) = (-l + v/TT2^)e-(-i+^'^)^V^-i+^'^)^dz. 

0<u<eq 

Plugging everything in, we can now say that 
Pxo{R{eq)>x) = 

+ 



g- 1 + VT+2i)3no g- (H-VT+2^)x 



g-(-l + /T+2^)3;og-(l-l-/T+2^)x 



1 



(-1 + VTT2^) 



and so after taking derivatives and multiplying by (—1), we find that the transient density of R{eq), 
for X > xq, is just 



fH(.,){x;x,) = (l + yTT2^)e-(-i+v™)-"e-(i+^)- 
1 



We will now focus on computing f R{e^)(x] xq) , for x < xq. After applying our weak- convergence 
results, we see that 



Pxo{R{eq) > x) 



]7r(dy) 



r ,.00 r , . .dPxni inf R(u) < z). 
Ey[e-i^»]TT{dy) °^o<u<e, ' - ' 



Evaluating this quantity, then taking derivatives shows that the transient density of R{eq) is just 

+ (1 + yT~r2^)e"("^+^^^)^"e"(^+^^^^^'' 



g - ( - 1 + X g - ( 1 + VT+2?) 2: 



(-1 + yTT2^)(2 - yTT2^)e-(-i+^"^)^«e-(i-^'^)^ 
y/l + 2q+l + q 



g- (- l+VT+2^)xo g- (1 + VT+2^)2; ^ 
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A Palm measures 



Throughout this paper, we assume that all of our random elements reside on a probability space 
(fi, J", P), where ft represents a complete, separable metric space, J" the Borel tr-field generated by 
the open sets of the metric, and P a probability measure on J^. These additional restrictions will 
be needed in order to properly define a collection of Palm measures, which are used to derive our 
main result. The reader should not be alarmed by such restrictions, as the space D[0, oo) endowed 
with the proper choice of Skorohod metric is a complete, separable metric space, and many queueing 
processes (and stochastic processes in general) can reside on such a space. Moreover, R+ is used to 
represent the nonnegative real line, and B the Borel tr-ficld generated by the open sets of M+. 

Let N := {N{t); t > 0} represent a point process on the nonnegative real line, with mean measure 
H, where fi{A) = E[N{A)] < oo for all bounded A G B. Under such assumptions, it is known that N 
induces a n-a.e. unique probability kernel V : K+ x F ^ [Oj 1]; where for each fixed E € T, VaiE) 
is a Borel measurable function in s, and for each fixed s € K.+ , Vs is a probability measure on 
The probability distributions of this kernel are referred to as the Palm measures of A'', and these are 
defined to be the measures that satisfy the following condition: for each B E B, and each A E J^, 



JB 

An important consequence of equation © is the Campbell-Mecke formula; see for instance Kallen- 
berg [25]. The proof of this formula follows from applying a monotone class argument to ([9]). 

Theorem A.l (Campbell-Mecke formula) For any measurable stochastic process {X[t)]t > 0}, we 
find that 



where Eg represents expectation, under the probability measure Vs- 

Throughout, we say that a stochastic process is measurable if it is measurable with respect to the 
(T-ficld A, which is generated by sets of the form A x C, where A G B, and C G F, i.e. if for each 



The Campbell-Mecke formula is a very important, fundamental result in the theory of Palm 
measures, and is typically the main tool used when applying Palm measures to a given problem. 
Readers wishing to consult a rigorous treatment of such measures are referred to Chapters 10-12 of 
[25] : other classical references on point process theory include the series of textbooks by Daley and 
Vere- Jones ^EM- 

A collection of sub-cr-fields {Fg; s > 0} of is said to be a filtration, if for each s < t, Fg C Ft- 
We say that a stochastic process {X{t);t > 0} is adapted to the filtration if, for each t >0, X{t) is 
measurable with respect to Ft- Associated with a filtration is a collection of cr-fields {Fs-; s > 0}, 
where Fs- is the smallest cr-ficld containing all cr-fields J>, for r < s. These are standard concepts 
within stochastic calculus, and can be found in virtually any textbook on the subject. Some examples 
of textbooks that focus on point processes, and include such concepts, are Bremaud [TT] and Baccelli 
and Bremaud [7]. 

We are now ready to quote a result that is used to derive the main result of this paper. Suppose 
N ■-= {N{t);t > 0} represents a point process on [0,00), and suppose {Ft;t > 0} represents a 
filtration, to which N is adapted. Within this framework, we say that N is an J^t-Poisson process, if 
(i) N is adapted to the filtration, and (ii) the distribution of N{a,b], conditional on Fa, is Poisson 
with rate 



for some deterministic function A : [0,cxj) — > [0, c») (i.e. N{a,b] is independent of Fa)- Under these 
conditions, we can apply the following result, which is a corollary of a time-dependent analogue of 
Papangelou's lemma for point processes; see [18] for details. 




(9) 



E 




BeB, {{t,uj);X{t,uj) eB}eA- 
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Proposition A.l If N is an Tt-Poisson process, then Vt = P on J-t-, for almost all t (w.r.t. 

Lehesgue measure). 
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